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Abstract 

Certain systems with slow driving and avalanche-like dissipation events are 

naturally close to a critical point when the ratio of two energy scales is large. 

The first energy scale is the threshold above which an avalanche is triggered, 

the second scale is the threshold above which a site is affected by an avalanche. 

I present results of computer simulations, and a mean-field theory. 
PACS numbers: 05.40.+j, 05.70.Jk, 05.70.Ln. 
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I. INTRODUCTION 



During the past years, systems which exhibit self-organized criticahty (SOC) have at- 
tracted much attention, since they might explain part of the abundance of fractal structures 
in nature Their common features are slow driving or energy input (e.g. dropping of 
sand grains |T]], increase of strain 0], tree growth 0, spontaneous mutations |Q) and rare 
dissipation events which are instantaneous on the time scale of driving (e.g. sand avalanches, 
earthquakes, fires, or a series of rapid mutations). In the stationary state, the size distribu- 
tion of dissipation events obeys a power law, irrespective of initial conditions and without 
the need to fine-tune parameters. There is, however, no reason to expect that systems with 
slow driving and instantaneous avalanches are always SOC. Such systems might also have 
many small avalanches which release only little energy, or only large avalanches which re- 
lease a finite part of the system's energy, or some combination of both. SOC systems are 
naturally at the critical point, due to a conservation law (sandpile model), a second time 
scale separation (forest-fire model), a competition between open boundary conditions and 
the tendency of neighboring sites to synchronize (earthquake model p|,||), or due to the slow 
driving alone ( "evolution" model). Often, the critical behavior breaks down when details 
of the model rules are changed (e.g. the boundary conditions in the earthquake model 
or the tree growth rule in the forest-fire model [^). 

There are certain parallels between these models and equilibrium critical systems, since 
both consist of many small units which interact with their neighbors, and since spin clusters 
in an Ising model or clusters of occupied sites in percolation theory can be compared to 
avalanches. However, the critical behavior of nonequilibrium systems can depend on mi- 
croscopic details, as mentioned above, in contrast to equilibrium critical phenomena, which 
commonly show universal behavior. Also, nonequilibrium systems do not satisfy a detailed- 
balance condition and can e.g. show periodic behavior. Furthermore, avalanches are usually 
released when some variable reaches locally a threshold, while other regions of the system 
might be far below the threshold, and consequently not all parts of the system look equal. 
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This can in particular result in more than one diverging length scale, as in the earthquake 
model or in the forest-fire model By contrast, in equilibrium systems the energy is 
an extensive variable, which means that all regions (which are large compared to the lattice 
size) are equal. This is e.g. the basis for hyperscaling relations. 

In this paper, I will add to the rich scenario of SOC systems by discussing models with 
slow driving and avalanche-like dynamics that show SOC behavior of a new kind. They have 
two threshold energies Ec and Eo and become critical with a power-law size distribution of 
avalanches in the limit of a separation of energy scales, {Ec — Eq)/ Eq oo. The size of the 
largest avalanche (i.e. the cutoff in the power law) depends on the degree of energy scale 
separation. These conclusions are derived from computer simulations in d = 2 dimensions, 
from intuitive arguments, and from analytical mean-field calculations. I will also argue that 
the critical behavior breaks down when certain details of the rules are changed, leading to 
synchronized dynamics. 

The outline of this paper is as follows: In section ||, the basic model is defined, and 
simulation results in two dimensions are presented. Section |T| defines related models that 
have the same behavior as shown by simulation results. Section |V| discusses the mean-field 
theory, and section |V| contains the conclusions. 

II. THE MODEL 

The model presented in this section is inspired by a deterministic forest-fire model in- 
troduced in a few years ago as a toy model for turbulence. Since that paper discusses the 
model only for one choice of the parameters, the generalization given below is natural and 
provides a deeper understanding of the model. The model is defined on a (i-dimensional 
hypercubic lattice with L'^ sites. In the beginning, each site is assigned a variable E which 
is randomly chosen from the interval [0, Ec) and can be interpreted as "tree height" or "en- 
ergy" . Sites with a height between and Eq are called "empty" , sites with a height between 
Eq and Ec are called "tree", and sites with a height above Ec are called "burning". All 
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heights are increased globally and very slowly. After some time, one of them reaches the 
threshold and becomes a burning tree. Fire spreads to trees on the 2d neighboring sites, 
and burning trees loose part of their height. After a fire has died, all heights increase until 
another tree reaches the threshold height, releasing a new fire. Driving is so slow that only 
a negligeable amount of energy is put into the system while a fire is spreading. In the sim- 
ulations, the system is not driven at all during a fire, and at the end of a fire the energy of 
all sites is increased by the difference between Ec and the largest height in the system. The 
precise rules for the spreading of the fire are the following: 

E — > E + Ec — Eq, if a nearest neighbor is burning and E e [Eq, Ec) . (2.1) 

E^{E-Ec)^r^: iiE>Ec. (2.2) 

They are apphed to each site at each time step, as long as there exist burning sites. The 
net effect of a fire on a tree of height E e [Eq, Ec) is a reduction in height to the new 
value {E — Eq)Ec/ {Ec — Eq). The interval [Eq, Ec) is linearly mapped on the interval [0, Ec). 
The distance in energy between neighboring sites therefore increases when both sites are 
burnt, making any kind of synchronization impossible for this model, and leading to chaotic 
behavior. A site can burn several times during a fire, until its height or the height of all its 
nearest neighbors is below Eq. We define the size s of a fire as the number of sites burnt 
during a fire, including multiple burnings. 

After a certain time, the systems reaches a stationary state with a certain mean energy 
and a certain size distribution of fires (averaged over either an ensemble of systems or many 
time steps). In this paper, we are only interested in the properties of this stationary state, 
and not in the transient behavior of the first time steps. Since the model is deterministic, 
chaotic, and dissipative, the stationary state is a strange attractor which occupies a vanishing 
fraction of the L'^-dimensional phase space (however, due to rounding errors in the computer, 
the model is not completely deterministic). Since the dynamics of the system are invariant 
when all energies (including the thresholds) are multiplied by the same factor, variation of 
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Ec alone gives all information about the system, and we can limit the discussion to the case 
Eq = 1. In the original version of the model, Eq = 1 and E^. = 2 0, and the size distribution 



of fires is a power law over a few orders of magnitude . 

When Ec is equal to i?o = 1, there are only fires of size one, since all neighbors of a site 
reaching the threshold are still empty. When E^ is slightly above Eq, sometimes neighbors 
of a tree reaching the threshold are between 1 and Ec, and fires burning several sites can 
occur. There is a cutoff in fire size which depends on Ec. With increasing Ec, the ratio 
between the number of trees and empty sites becomes larger, and the cutoff in fire size 
increases. We expect therefore fires spreading over the whole system when Ec is larger than 
some critical value E*. The simulation results, however, do not show such a transition for 
any value of Ec that can be studied with linear system sizes L < 400. Fig. [l| shows the size 
distribution n{s) of fires for L = 400. The simulations were started with a random energy 
distribution. To assure that the system is in a stationary state, the first 10000 simulation 
steps were discarded, and n{s) was determined from the following 10000 steps. The fire 
size distribution decays as with r ~ 1.07 ± 0.03, and with a cutoff for larger fire sizes. 
Simulations for Ec > 3.0 show finite-size effects since they lead to avalanches of the order of 
the system size L^. This explains the bump in the fire size distribution for Ec = 4.0. Similiar 
bumps occur already for smaller values of Ec when the system size is smaller. Of course, such 
bumps should also occur when the system is beyond a critical point and has true infinite 
avalanches. But the simulations do not show the other signatures expected for a phase 
transition: Beyond a critical point, the bump size increases with increasing distance from 
the critical point, and the cutoff in cluster size diverges when E* is approached from below. 
Although it is impossible to derive from our simulation results the analytical dependence of 
the cutoff cluster size on Ec, we can clearly see that it has no tendency to diverge when Ec 
increases. Although these simulation results allow no conclusive prediction whether there 
exists a finite E* that is much larger than the values studied in the simulations, it seems well 
possible that the critical energy is infinity. This hypothesis is supported by the probability 
distribution of energy in the system, shown in Fig. ^. This probability distribution has a 



rapid decay for energies larger than Eq = 1, and seems to approach a constant form as Ec 
increases. This means that the fraction of excitable sites is almost the same for all energies, 
indicating that there is no phase transition in the system. A simple consideration shows 
that the mean fire size s = sn{s) diverges in the limit —> oo if the function P{E) 
approaches a constant form: Let v be the (very slow) driving velocity. Per unit time dt, 
L'^sP{Ec)vdt sites participate in a fire. If the function P{E) approaches a fixed function 
in the limit Ec — > oo, this number must also be independent of Ec. Consequently, sP{Ec) 
must approach a constant for Ec oo. If P{E) changes only slightly with increasing Ec, 
only a vanishing fraction of sites can have large energies, and P{Ec) for Ec oo. This 
implies that s oo. 

The following intuitive argument gives a deeper understanding why E* cannot be finite 
in our model: Assume that there is a finite E*, above which the mean fire size becomes 
proportional to the system size, s oc L^. These infinite fires dissipate a finite fraction of 
the system's energy. Consequently, a finite fraction of the system's energy must be put into 
the system between fires. This is only possible if the highest energy in the system after an 
fire is at a finite distance from Ec. However, the rules of our model do not allow such an 
energy gap, since the energy of sites that are very close to Ec is decreased only by a very 
small amout if this site is involved in a fire. Since this argument is very general, we expect 
a similiar behavior also in higher dimensions. In the language of critical phenomena this 
means that the lower critical dimension of the model is infinity. 

The scenario in one dimension is slightly different. Since any site with an energy below 
Eq = 1 stops the fire, less and less sites have energies below Eq as Ec increases and the 
mean fire size becomes larger. The function P{E) therefore does not approach any finite 
limit function as Ec diverges. 

Once we have understood the mechanism leading to -E* = oo, we can easily construct 
other models with the same behavior. The next section introduces two of these models 
and shows simulation results in two dimensions. Section |V| then presents the mean-field 
theory, showing explicitely that E* = oo. Since the mean-field theory usually describes 
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well systems in high dimensions, this analytical calculation supports the conclusion that the 
behavior observed in two-dimensional simulations persists in higher dimensions. 

III. RELATED MODELS 

The above model is difficult to investigate analytically, since the energy of a site after 
a fire is related to the energy before the fire in a nontrivial way, leading to complicated 
correlations. Models where the energy after a fire is essentially determined by a random- 
number generator are easier to understand. Let us consider therefore first a model that has 
the following version of rule 2 (Eq. (|2.2|) ) above: 

E — > EcX rand , ii E > E^ . (3.1) 

rand is a random number equally distributed in the interval [0, 1). Like the previous model, 
this model does not allow an energy gap between the highest energy and E^., and we expect 
again E* = oo. Since the fire jumps back and forth between neighboring sites until at 
least one of them is below Eq = 1, a site that has participated in a fire is below Eq with a 
probability larger than 0.5 even for large Ec, suggesting again a fixed limit limE,.^oo P{E), 
with only very few sites with high energies. Since all energies between zero and Eq are 
chosen with equal probability according to rule ( |3.1|) , we expect that P{E) is linear in the 
interval [0,i?o], as can also be seen from the explicit calculations in the next section. I 
performed the simulations under similar conditions as for the previous model. The resulting 
distributions n{s) and P{E) are shown in figures ^ and ^. The exponent r characterizing the 
fire size distribution is r = 1.08 ± 0.02. This value is compatible with the one obtained for 
the deterministic model, and the two models possibly belong to the same universality class. 
The probability distribution of energies P{E) is linear for E < Eq (this cannot easily be seen 
in this logarithmic plot), and seems then to dacay exponentially, reaching a constant value 
proportional to 1/{Ec — Eq), a behavior that will also be obtained from the mean-field theory 
below. In order to make the interesting part of the curves more visible, the horizontal part 
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has been cut off for energies larger than 3. Fig. |^ shows the mean fire size s as function of Ec- 
The best fit is obtained from a stretched exponential s oc exp[c{Ec — Eq)'^], with c = 2.4 and 
7 = 0.37. We will see below that the mean-field theory gives a simple exponential increase. 
Fig. ^ shows a snapshot of part of the system for Ec = 5.0. All sites with E < 1 are black, 
other sites have a grey shade that indicates the energy. The scale for the grey shade is not 
linear. 

A model that is closer to the deterministic model of the previous section and that can 
also be analyzed in mean-field theory has the following rule 2: 

E — >{E -Ec + Eq) X rand , ii E > E^ . (3.2) 

With this rule, the energy of a site is always decreased during a fire. Nevertheless, the 
probability that a given site will have an energy close to E^ after a fire does not vanish, and 
we expect again E* = oo. The fire size distribution for this model is shown in figure ^. The 
cutoff fire size increases much faster than in the previous model, and we cannot perform 
simulations with Ec > 2.5 without getting finite-size effects. This is because large energies 
occur now with a much smaller probability than before, reducing P{E^. If we assume that 
the fraction of sites below the threshold Eq = 1 of the same order in different models, 
it follows immediately that s (which is proportional to l/P{Ec) as we have seen before) 
increases. Unfortunately, we cannot even extract a useful information about r from the 
simulation results, since the system is still far from the asymptotic behavior. 

IV. MEAN FIELD THEORY 

The two models defined in the previous section can be solved analytically in mean-field 
theory. In a mean-field theory, all correlations between different sites are neglected, and the 
density P{E) is calculated self-consistently. A random-neighbor version of the above models 
shows the same behavior as the mean-field theory. In this version, the z = 2d neighbors of 
a burning site are chosen at random for each simulation step. In the original models, the 
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fire jumps back and forth between neigliboring sites, until at least one of the two has an 
energy below Eq. We take this into account in the mean-field theory by requiring that a 
finite fraction a of all sites that participate in a fire end up with an energy below Eq, where 
a is independent of Ec. The mean-field equations for the first stochastic model that assigns 
a random energy from the interval [0, E^) to each burning tree are the following 

^ - + vasPiE.yEQ iorE<EQ.l; 

dP{E,t) ^ dP{E,t) vsPjEct) vP{E,t)P{Ec,t){s-l) 

dt Ec-Eq jg^P{E',t)dE' 

ioiEQ<E<Ec. (4.1) 

The first term on the right-hand side is due to the driving of the system that replaces P{E, t) 
by P{E — f dt) after time dt. The second and third term describe the change in P{E, t) due 
to fires. The fraction of trees that catch fire per time dt is vdtsP{Ec). The fraction a of 
the burning sites become evenly distributed in the interval [Oj-Eq]? while the rest is evenly 
distributed in the interval [Eq, Ec). The fraction of sites with E G {Eq, Ec) that catch fire per 
time dt is given by the third term in the second equation. Since there are no correlations in 
the system, the probability that a site of energy E becomes involved in a fire is proportional 
to the number of sites with this energy. 

Since we are interested in the stationary state, we can drop the time-dependence and set 
dP/dt = 0. The solution of Eqs (^4.1| ) is then given by 

P{E) = aP{Ec)sE/EQ for E < Eq- 

exp [-P{Ec){-s-1){E-Eq)II3] 
for Ec> E>Eq, (4.2) 

with p = /I; P{E')dE'. Normalization requires P = 1- aP{Ec)EQs/2. We still have to 
find P{Ec) and s to complete the solution. Since we are interested in the critical behavior 
of the model, let us focus on the situation where s becomes very large. The mean fire size 
s diverges as [{l/z) — P]^^ when /? approaches the percolation threshold 1/z. This allows 
us to replace [s — 1) with s, and the right-hand side of Eqs (f4.2|) depends only on the 
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product P{Ec)s, but not on these two quantities separately. Since P{E)dE approaches 
a constant {1 — 1/z) as s diverges, this product must become also a constant C for large s. 
From (3 = 1 — CaE^jl = 1/z follows C = 2{z — 1)/ zaE^. The solution of our mean-field 
equation therefore becomes (for large s) 



P{E) = CaE/Eo for E<Eo; 

:i-a) 



Ca- 



exp[-Cz{E - Eo)] for E>Eo. (4.3) 



ziE, - Eo) 

This solution shows that P{Ec) vanishes in the limit Ec — > oo, and we conclude that s 
diverges exponentially as exp[Cz{Ec — Eo)]. The exponent 7 charecterizing the stretched 
exponential increase of s is 7 = 1 in mean-field theory. The exponent r characterizing the 
size-distribution of fires is r = 1.5 in mean-field theory, since the density of tres with E > Eo 
is at the (mean-field) percolation threshold. (For the mean-field theory of percolation see 
e.g. Q.) 

The mean-field theory for the second model is slightly different. While the solution for 
E < Eq has the same form as before, the probability distribution for energies larger than 
one is obtained from the equation 

d-P(-g) _ ^ . sP{E,) P{E')dE' P{E)P{E:){s - 1) 
AE ' Jg^ J^r PiE")dE"dE' !gP{E')dE' 

Introducing R{E) = J^" P{E')dE', this becomes 

d^R dR ,^ , 
— +a-+bR^O, (4.5) 

with a = P{E^){s - l)/R{Eo) and b = [1 - a)sP{E^)/ jg' R{E)dE. Solving this equation, 
and calculating P{E) = —dR/dE gives 

P{E) = Aexp[-Ai(E - Eo)] + B exp[-X2{E - Eo)] , (4.6) 



with A1.2 = (a/2) ± ^ (ct^/4) — b. The coefficients A and B have to be determined by 
the condition that P{E) is continuous at E = Eo, and by the normalization condition 
J^" P{E)dE = 1. As before, the P{E) depends only on the product P{Ec)s for large s. 
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which must become a constant C — 2{z — l)/azEo as s diverges. Again, we focus on the 
hmit of large s, where a — Cz. We want to show that s can only diverge in the limit 
Ec ^ oo, but not for a finite Ec = E*. For this purpose, let us assume a finite E* and show 
that this leads to contradictions. The conditions P{Ec) — ai, Ec — E^ and P{Eq) — Ca 
lead to 



sinhUEc- E)J(ay4) -b] 

P(E) = CaeM-a{E - ^o)/2] — for Ec = E: . 

smhiiEc - Eo)^{ayA) - b] 

It is easy to see that the normalization condition is J^^ P{E)dE = l/z cannot be satisfied. 
We have 



P(E)dE < aC/a ^a/z <l/z, 



since a < 1. 

Consequently, s can only diverge in the hmit Ec — > oo. In this limit, J^^ R(E)dE 
diverges as fast as Ec — Eq since most sites have energies much smaller than Ec- From 
the normalization and continuity of P{E) we find then to leading order in l/{Ec — Eq): 
A ^ aC - B, B ^ {1 - a)/zQ{Ec - Eq), Xi ^ Cz - A2, and A2 = l/Q{Ec - Eq). The 
constant Q depends on a and is defined by the equation 

Sg + (1 - Q) exp(-g-i) = 1/(1 - a) . 

As in the first model, the asymptotic form of P{E) is a simple exponential decay for E > Eq 
with a decay constant Ai = 2(2; — Vj/aE^. 



V. CONCLUSION 

In this paper, I have shown that certain systems with avalanche-like dynamics and two 
enery scales have a power-law size distribution of avalanches over many orders of magnitude 
when these scales are widely separated. For some of these models, a very moderate scale 
separation {Ec — Eq)/Eq ~ 3 gives already scaling over more than four decades. Since a 
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separation of enery scales can occur naturally, it can give rise to power laws in nature. These 
models belong therefore to the class of SOC systems. They are to some degree similiar to 
the SOC forest-fire model. In that model, a separation of time scales for tree growth and 
lightning strokes is responsible for the observed power laws. When the lightning probability 
is very small compared to the tree growth probability, an isolated tree has to wait for a long 
time before it catches fire spontaneously. The separation of energy scales in this paper has 
a similiar effect: When {Ec — Eq)/Eo is large, a newly grown tree (energy Eq) that is not 
ignited by its neighbors has to wait for a long time until it reaches the threshold Ec and 
catches fire. But there are also major differences between the two models. The models of 
this paper have a continuous energy scale, and a finite fraction of all burning trees remain 
trees after a fire. The exponent r characterizing the size distribution of fires seems to be 
different. 

Therefore, the models studied in this paper are a new class of SOC systems. Many 
variations of these models are expected to be also SOC, as long as trees with an energy close 
to Ec have a nonvanishing probability of loosing only a small fraction of their energy during 
a fire. However, if burning trees turn always to empty sites, the system becomes completely 
synchronized when Ec > 2Eo. SOC systems seem to be in some cases systems that fail to 
synchronize, as also observed in |^ for the earthquake model. 

Let us conclude with a caveat: Not all models that satisfy the above-mentioned critieria 
are SOC. To illustrate this, let us consider a model where a burning tree reduces its energy 
by some small amount AE ^ Eq with some small probability g, and becomes an empty site 
with zero energy with probability 1 — g. If we begin with a completely synchronized system, 
after the first fire a fraction 1 — g of all sites will have zero energy, and the remaining sites 
are very close to the threshold. All these sites will burn down before the other sites reach 
the energy Eq, and there will again be a fire that spreads over the entire system. Thus, the 
conclusion of this paper illustrates what is already said in the introduction: Slowly driven 
non-equilibrium systems with avalanche-like dynamics show a richness and complexity of 
behavior that still has to be fully explored. 
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FIG. 1. Size distribution of fires for L = 400 in the deterministic model. 
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FIG. 2. Probability distribution of energies for L = 400 in the deterministic model. 



15 




s 



FIG. 3. Size distribution of fires for L = 400 in the first stochastic model. 




FIG. 4. Probabihty distribution of energies for L = 400 in the first stochastic model. 
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FIG. 5. Mean fire size s as fuction of the critical energy Ec- The soUd fine is a stretched 
exponential s oc exp[— 2.4(£^c — 1)*^'^^]. 




FIG. 6. Shapshot of the first stochastic model. Sites with E < Eq are black; the energy of the 
other sites is indicated by the grey shade. 
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